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The dc Josephson effect in a one-dimensional Tomonaga-Luttinger (TL) liquid is studied 
on the basis of two bosonized models. We first consider a TL liquid sandwiched between two 
superconductors with a strong barrier at each interface. Both the interfaces are assumed to be 
perfect if the barrier potential is absent. We next consider a TL liquid with open boundaries, 
weakly coupled with two superconductors. Without putting strong barriers, we instead assume 
that the coupling at each interface is described by a tunnel junction. We calculate the Josephson 
current in each model, and find that the two models yield same results. The Josephson current 
is suppressed by repulsive electron-electron interactions. It is shown that the suppression is 
characterized by only the correlation exponent for the charge degrees of freedom. This result 
is inconsistent with a previously reported result, where the spin degrees of freedom also affects 
the suppression. The reason of this inconsistency is discussed. 
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1. Introduction 

Progress in micro-fabrication technology has enabled 
us to prepare normal-conductor-superconductor (NS) 
composite systems without a Schottky barrier at the 
NS interface. 1 A highly controllable normal segment has 
been realized by using a two-dimensional electron sys- 
tem formed at semiconductor heterostructures. 2 We are 
interested in how the electronic transport properties are 
influenced by superconducting proximity effect. Partic- 
ular attention has been focused on the low-temperature 
transport properties with the phase coherence of elec- 
trons in the normal segment. Various experiments have 
been performed to reveal unusual features of the con- 
ductance in NS systems and the Josephson effect in SNS 
systems. 2 ' 3 

Progress in micro-fabrication technology has also made 
it possible to prepare a very narrow quantum wire which 
can be regarded as an idealistic one-dimensional (ID) 
electron system. 4,5 The electron transport in a ID sys- 
tem with a few barriers has been studied extensively. 6 
The central problem addressed there is how mutual 
electron-electron interactions affect the transport prop- 
erties. Electron-electron interactions greatly affect the 
low-energy properties of a ID electron system. A ID 
interacting electron system with a gapless excitation is 
generally described as a Tomonaga-Luttinger (TL) liq- 
uid. 15-17 The strength of electron-electron interactions is 
characterized by the correlation exponents K p and K a , 
where K p (K a ) is related to the charge (spin) degrees 
of freedom. For example, K p < 1 and K a — 1 cor- 
responds to the spin-independent repulsive interaction 
cases, and the system is reduced to the noninteracting 
electron gas when K p — K a — 1. Reflecting the nature 
of the TL liquid, we expect that the electron transport 
in a ID system shows an anomalous property, which is 
not seen in a conventional Fermi liquid. Kane and Fisher 
studied the electron transport in a spin-less TL liquid 
with a single barrier potential of 5-function type. 6 Using 
rcnormalization-group argument, they found that for re- 



pulsive electron-electron interactions with K p < I, the 
potential becomes large after renormalization even if it 
initially is very small. This means that electrons at low 
energy see the potential as if it were effectively infinite. 
Thus, the low-energy properties are described by the 
open-boundary fixed point, at which the ID system is ef- 
fectively disconnected at the barrier. This anomalous be- 
havior strongly affects the transport through the barrier. 
It is shown that the conductance G is suppressed with 
decrease of temperature T as G cx T 2K p -2 , and eventu- 
ally vanishes at T = 0. This ia a clear manifestation of 
the TL-liquid behavior. The extention to spin-dependent 
cases is achieved by Kane and Fisher, 8 and Furusaki and 
Nagaosa. 10 The conductance in this case also obeys a 
power-law behavior as a function of T, and its exponent 
is determined by K p and K a , Furthermore, the electron 
transport in a TL liquid with a double-barrier structure 
was also studied by these authors. 7-9 

Inspired by the developments mentioned above, con- 
siderable attention has been attracted to the transport 
properties of a TL liquid coupled with superconduc- 
tors. 18-27 The central problem is how the supercon- 
ducting proximity effect in a ID system is modified 
by electron-electron interactions. Particular interest is 
focused on the Josephson effect in a superconductor- 
TL-liquid-superconductor (STLLS) system, 19-21 ' 24 - 26 ' 27 
which is the issue of the present paper. 

Let us consider a clean ID noninteracting electron 
system of length L sandwiched between superconduc- 
tors (see Fig. 1(a)). We focus our attention on the long- 
junction case, where v-p/L (vp: Fermi velocity) is much 
smaller than the energy gap A of superconductors. In 
the low-temperature regime of T <C v-p/L, the Joseph- 
son critical current obeys j c oc evp/L. 2S In the case of 
an STLLS system, we expect that the critical current 
deviates from the inversely linear L-dependence due to 
the TL-liquid behavior. This problem is first discussed 
by Fazio et a/. 19 ' 20 They considered a clean TL liquid of 
infinite length weakly connected to two superconductors 



2 J. Phys. Soc. Jpn. 



Full Paper 



Yositakc Takane 



s l 



Sr 



X=0 t v X=L 

(a) 



x=0 x=L 

A — A 



(b) 



Sr 



Fig. 1. The geometries discussed in the text, (a) One-dimensional 
electron system of length L sandwiched between two supercon- 
ductors, (b) Onc-dimcnsional electron system of infinite length 
weakly coupled with two superconductors. Sl (Sr) represnts the 
left (right) superconductor. 



(see Fig. 1(b)), where the separation between the two 
NS contacts is L. They found that when T <C Vf/L, 
the L-dependence of the critical current is given by 
j™° oc (l/L) K p + K "- 1 . This power-law L-dependence 
is a clear manifestation of the TL-liquid behavior. It in- 
dicates that as long as K~ x + K a > 2, the critical cur- 
rent is suppressed by electron-electron interactions com- 
pared with the noninteracting case. This result shows 
that both the charge and spin degrees of freedoms in- 
fluence the suppression of the critical current. It should 
be noted that in calculating the critical current in the 
STLLS system, Fazio et al. completely neglected poten- 
tial scattering at the contacts between the TL liquid 
and the superconductors. Fabrizio and Gogolin 29 pointed 
out that since the potential induced by the contacts be- 
comes large after renormalization even if it initially is 
very small, 6 the TL liquid is finally disconnected at the 
NS contacts. They suggested that in such a situation, we 
should employ a TL liquid with open boundaries to cor- 
rectly take account of the potential scattering in studying 
the Josephson effect. Maslov et al. 21 studied a different 
model, in which a TL liquid of length L is sandwiched 
between superconductors (see Fig. 1(a)). They presented 
a uscfull bosonized description of the STLLS system, and 
calculated the Josephson current when the NS interfaces 
are perfect. They also examined the L-dependence of the 
critical current in the presence of a strong barrier at each 
NS interface. However, they do not correctly take ac- 
count of NS boundary condition in their renormalization- 
group argument. The works mentioned above treat a TL 
liquid with repulsive electron-electron interactions. The 
Josephson effect in attractive interaction cases has been 
studied extensively by Affleck et al. 2e 

In this paper, we study the dc Josephson effect in a 
TL liquid using two different models. We first consider 
a TL liquid of length L sandwiched between two super- 
conductors with a barrier at each NS interface. Both the 
NS interfaces are assumed to be perfect if the barrier po- 
tential is absent. We call it the interface-barrier model. 



Particular attention is focused on the limit where the 
barrier potential is very strong. This model is equivalent 
to that treated by Maslov et al. 21 We next consider a TL 
liquid of length L with open boundaries, which is weakly 
connected at its left (right) end with the left (right) su- 
perconductor through a tunnel junction. We call it the 
weak-coupling model. This model is equivalent to that 
suggested by Fabrizio and Gogolin. 29 It is shown that 
the two models provide us essentially the same expression 
of the Josephson current in both the high-temperature 
regime of Vf/L <C T <C A and the low-temperature 
regime of T < v F /L. We find that j c oc (l/L) 2 ^ 1 " 1 
in the low-temperature regime. Clearly, this result is in- 
consistent with Fazio et al.'s result, 19,20 where the in- 
dependence is characterized by both K p and K a . The 
incosistency is attributed to the difference in the treat- 
ment of the potential induced by the NS contacts. The 
basic assumption of Fazio et al. is that the supercon- 
ductors do not influence the uniformity of the poten- 
tial along the ID system, while we assume that the low- 
energy property of an STLLS system is described by the 
open-boundary fixed point. Since the potential induced 
by the contacts becomes large after renormalization even 
if it initially is very small, Fazio et a/.'s result is unstable 
against the renormalization unless the induced potential 
is completely negligible. If the induced potential cannot 
be neglected, the ID system flows to the open-boundary 
fixed point. That is, the ID system is effectively discon- 
nected at the NS contacts. Fazio et al.'s model is not 
appropriate in this situation, and we expect that our 
weak-coupling model provides us correct results. 

This paper is organized as follows. In the next sec- 
tion, we present the interface-barrier model on the ba- 
sis of a bosonized description of an STLLS system. The 
bosonization procedure based on Haldanc's argument is 
outlined in Appendix A. In §3, we calculate the Joseph- 
son current in the interface-barrier model using the in- 
stanton approximation. The Josephson current in the 
weak-coupling model is studied in §4. We employ the 
bosonized description of a TL liquid with open bound- 
aries given by Fabrizio and Gogolin, 29 and calculate the 
Josephson current within the lowest-order perturbation 
with respect to the coupling between the TL liquid and 
the superconductors. Section 5 is devoted to discussion. 
The argument concerning to the interface-barrier model 
has been briefly reported in ref. 24. Some errors in ref. 24 
are corrected in this paper. We set h = ks = 1 through- 
out the paper. 

2. Interface-Barrier Model 

We consider a ID electron system of length L sand- 
wiched between superconductors (see Fig. 1(a)). The su- 
perconducting order is characterized by the pair poten- 
tial A(x). We assume that A(x) is given by 28 

( Ae ixi for x < 
A(x) = I for < x < L (1) 
[ Ae iX2 for L < x, 

where xi (X2) is the macroscopic phase of the left (right) 
superconsuctor. We assume that the length L is much 
longer than the coherence length which is defined by 
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£ = vp/A with the Fermi velocity v F . A quasiparticle 
state in our system is described by the Bogoliubov-de 
Gennes (BdG) equation. 30 Due to Andreev reflection at 
NS interfaces, quasiparticle states depend on the phase 
diffcrncc \ = X2 — Xi- We introduce a set of eigenfunc- 
tions of the BdG equation. The eigenfunctions are easily 
obtained if we assume that perfect Andreev reflection 
(i.e., zero normal scattering) is achieved at the NS inter- 
faces. In the perfect Andreev reflection case, we obtain 
the low-energy eigenfunctions in the Nambu notation 

e -i?i +e ig + (iE+?/2) 



Ug(x) 


\ 1 












/ V2L 



Ak-px 



(2) 



e iil- e -iq-(x+^/2) 



-ik-px 



where q = ir(n 
and 



l/2)/(L + with n = 0,±1,±2, 



q± = q± 



x 



2(i + 



Xi 

4 + 2 ' 



(3) 

(4) 
(5) 



We have assumed that the Fermi wave number k F is given 
by fcp = irno/{L + £), where no is an integer. Let e q+ and 
e q - be the linearized eigenenergy of the state eq. (2) and 
that of the state eq. (3), respectively. They are given by 

(6) 



e q ± = v F q ± 



1 

2m 



X 



In deriving eqs. (2) and (3), we have used the following 
approximation 

e q ± + i^A 2 -e 



4± 



expi 



f- 

V2 



A J ' 



(7) 



which is justified when e q ± -C A. 

We assume that any quasiparticle state in our system 
is described by eqs. (2) or (3). This is the key approxima- 
tion allowing us to perform the bosonization procedure. 
Clearly, this cannot be justified for q of the order of, or 
greater than, However, high-energy states compared 
with A play only a minor role in low-energy properties, 
which are of interest to us. To study Joscphson effect, 
we must take account of the zero modes 17 as well as 
the nonzero modes. To do so, we introduce the winding- 
number operators J and M, and the zero-mode opera- 
tors ip p and <p a , which satisfy [J, ip p ] — [M, tp„] = 2i and 
[J, if a ] — [M, <fp] = 0. The operator J (M) is related to 
the excess charge (spin) accumulated in the ID system. 
The nonzero modes are described by the boson opera- 
tors a q , a q , (3 q and /?|. The bosonization procedure is 
outlined in Appendix A. 

In terms of the operators introduced above, the 
bosonized Hamiltonian is given by 



IE (* 



+ X] { V P<l a \ a 1 + V <y<l h \ h q) . ( 8 ) 
q>0 

where v p (v a ) is the velocity of charge (spin) excitation 
and q = n(n+l/2)/(L+£) with n = 0, 1, 2, • • • . In eq. (8), 



K p and K a are the correlation exponents, which charac- 
terize the interaction strengths. For example, K p < 1 
and K a = 1 corresponds to the spin-independent repul- 
sive interaction case, and the system is reduced to the 
noninteracting electron gas when K p = K a = 1. It is 
important to note that eigenvalues of J + M are limited 
to be even integers, as shown in Appendix A. 

We express the electron- field operator as tp s (x) , where 
s = ± denotes spin index, and decompose it as 



1p s (x) = 1pu(x) + 1p2s(x), 



(9) 



where Vis and \p2s represent the right and left movers, 
respectively. To express Vis and V2s j we define the phase 
fields: 31 



e + (x) = ipp + 9 + (x), 



(10) 

(ii) 



jm(x + L+^)^ <> (,r). (12) 



$_(x) = tp a +<j>_(x). 



(13) 



The nonzero-mode components 0± and </>±, which de- 
scribe low-energy fluctuations around the zero modes, 
arc expressed in terms of the boson operators as 



-(*) 



-(*) 




e aq ' 2 cos 




e aq/2 smq(x+i 



Q (* + 2 ) ^ ~ ^ ' 
(14) 

(4 + a i) ' 



P q>0 



<j> + {x) 



cj>_{x) = 





e aq/2 sin q I x 



e ' aq/2 cosqlx 



(15) 



(16) 



q>0 



I) (b\-b q ), 



(17) 

where a is a positive infinitesimal. The phase field 6+ is 
regarded as the phase of the charge density wave, while 
6_ corresponds to the Josephson phase. The phase fields 
$ + and $_ are related to the spin degrees of freedom. 
Using the phase fields, we can express Vis and ip 2s as 32 

ipi s (x) = 1. — cxp (ik F x + ^(0+(x) +Q-(x) 

+ s$+(x) + «$_(a0)^, (18) 
exp I - ik F x + ^ (-9+ (x) + 6_ (x) 



lp2s(x) = 



V2 



- s$ + (x) + s$-(x)) 



(19) 



Now we put a barrier at each NS interface. For sim- 
plicity, we assume that both the barriers have the same 
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strength. Then, the corresponding Hamiltonian is 

# B = -Vu W(^(°hM0) + h.c.) 

+ (^(L)fe(I)+h.c.)j. (20) 

To express H B , it is convenient to introduce the new 
phase variables: 8,9 



= 1 (<?+(£) + MO)), 

9 = (9 + (L) - 9 + (0)) , 

4> =\{4>+{L) + <!>+$))■ 
4> = {<f>+{L) - <f>+(0)) . 



(21) 
(22) 
(23) 
(24) 



The time-derivative of 9 (4>) is related to the charge (spin) 
current, while 9 (0) is related to the excess charge (spin) 
accumulated in the ID system. In terms of them, we find 
that 



Wo 
ira 



cos (i^p + 9 + kpL) cos ( - + k-pL 



X cos 



+ sin (ip p + 9 + kpL) sin I - + kpL 



x sin ^</> + -ttM j sin + ttM 

'(25) 

where we approximated as (3irM/2)(l + £/(3L)) ~ 
3ttM/2. 

3. Josephson Current in the Interface-Barrier 
Model 

In this section, we consider the interface-barrier model 
in the limit where the barrier potential is very strong. 
The total Hamiltonian of our system is H tota \ = H + 
H B . We calculate the partition function Z in terms of an 
imaginary-time path integral. The dc Josephson current 
is obtained through the well-known relation 

, x d\nZ , x 

j( X ) = - 2eT —- ( 26 ) 

Since M commutes with H B , then M remains a con- 
served quantity despite the presence of the barriers at 
the NS interfaces. By contrast, J is not a good quantum 
number when V ^ 0. The partition function is expressed 
as 



M = 
(M=even) 



7 ( R™° M 

Z= E ^l^l^ 

M = — na x 

Vf = -oo \ ' 



where 

z ± = jr r ivp/n 





Va\Va q Vb\Vb q 



xcxpiiy dr^^drifip - S z - S NZ - S B 



with C+ = and £_ = w, and 

S z = [ 
Jo 



L 
dr- —(d T ip p ) 2 , 



4-KVpKp 



(28) 
(29) 



S NZ = (aj5 T a 9 + v p qa\a q 

Jo q>0 

+ b\d T b q + v a qb\b q ) 1 (30) 

rP 

S B = / drF B (r). (31) 
Jo 

In eq. (28), the integration over tp p is carried out un- 
der the boundary condition of (p p ((3) — f P (0) = 2-kvh. 
This boundary condition is imposed to take the discrete 
nature of J into account. Furthermore, we used (± to 
maintain the constraint J + M = even. 33 

Before considering the strong-barrier limit, which is of 
interest to us, we treat the case of Vo — ► to examine 
the validity of eq. (27). When S B = 0, _thc partition 
function factrizes as Z = Z(x)Z, where Z{\) and Z are 
the contribution from the zero modes and that from the 
nonzero modes, respectively. The former contribution Z 
is given by 

z= ex p -0 ■-- - I Z + 



M = — oo 
(M =cvcn) 



IK„L 



E ex P (-/3^]2 . ,32, 



(M = odd) 



where 

oo 

= E 



m— — oo 



x cxp 



dT*±^d T <p p -sA . (33) 



The latter contribution Z, which docs not depend on x, 
is irrelevant for the Josephson effect. We decompose ip p 
as (p p (r) = (p p o + 27rmr//3 + <p p (r) with <p p (j3) = <p p {0). 
Substitution of this into eq. (33) yields 



oo r 

z±= E e M-^hr m2 + Mx + C±))- (34) 



m— — oo 



Applying Poisson's summation theorem to eq. (34) and 
then substituting the resulting expression into eq. (32), 
we obtain the correct partition function in the absence 



M = - 
(M=odd) 



J. Phys. Soc. Jpn. 



Full Paper 



Yositake Takane 5 



of the barriers, 21 



J.M ^ 



J,M 
( J+M = even) 



(35) 



where the summation over J and M is carried out under 
the condition of J + M = even. 

Now we turn to the strong-barrier limit, where H B is 
regarded as a strong pinning potential for the phase fields 
9, 9, 4> and cp. It is convenient to introduce the effective 
action which satisfies 

J n g>0 Va\Va q Vb\Vb q exp (-S NZ - S B ) 

jn q>0 ValVa q VblVb q exp(-S™) 

J V6V6V4>V4>exp (-S^ - S B ) 



J V6V8V<j>T><t> exp (-S' N ff z ) 
This action is given by 



(36) 



UJ 

+ 8^bS^ 1 *)^M, (37) 



where 



Jp(w) = 



cosh ( ) + cosh 

V 



J p {w) 



2cosinh (^) 1 """" V v p ) \ v p J 
cosh + cosh 



2o> sinh 





(39) 



■7„(w) 



2o> sinh 



' 1 cosh (^)- cosh 



inh (-) 



cosh [ — 1+1 



(40) 



-7„M = 



1 



2a, sinh (^) . 

+ C ° Sh (S) 



ro S h(^Vc OS hf (L -^ 



(41) 



The derivation of the effective action is outlined in Ap- 
pendix B. From eqs. (40) and (41), we see that J" 1 ^) 



and J~ 1 (w) are of the order of A in the limit of uj — > 0. 
This indicates that (f> and <f> have a mass gap of the or- 
der of A and thereby their low-energy fluctuations are 
strongly suppressed. We thus see that 4> and <p are pinned 
by H B when T< A. This is an important feature of the 
phase variables in our problem. We can understand the 
feature by noting that only a spin-singlet pair of two 
electrons can transfer the NS interfaces due to Andreev 
reflection. Detailed discussion on the characteristic fea- 
tures of Sftf is given in §5. 

Since cj> and <fr are pinned by the barrier potential, we 
can set <fi = <fi = without loss of generality. Further- 
more, we note that M plays only a minor role in chang- 
ing the location of the potential minima, at which ip + 9 
and 9 are pinned. Consequently, S B can be simplified to 



S B = - 



4V 



TXOL 



Jo 



dr cos(v?p + 9) cos ■ 



(42) 



where k?L = (mod n) is assumed for simplicity. Since 
4> and 4> are irrelevant for our argument, S^g is reduced 
to 

1 £J-V)0>)0(-^) 



cNZ 
°cff 



2nK p (3 ■ 



1 



8nK p (3 



(43) 



It should be noted here that the high-frequency cutoff 
A/2 (2 A) must be introduced for the summation in the 
first (second) term to avoid ultraviolet divergence (A is of 
the order of A). To do so, we introduce an additional ac- 
tion which serves as the high-frequency cutoff for S^, 10 



'^M'+tK, («) 



where M w 2/A and M w 1/(2A). After these treat- 
ments, eq. (28) is reduced to 



Z± 



E 



3 im(x+C±) 



J Aip — 27rm J 



V9V9 



exp (-S* 2 



' A(p p — 2nm 

S™ z - S A - S B ) 



(45) 



In the strong-barrier limit, the change in ip p with m ^ 
occurs in collaboration with instanton tunneling of 9 
and 9 from a potential minimum (ip p + 9,9) = (kir, 2lw) 
to an adjacent minimum {tp p + 9,9) = ((k ± l)n, (21 ± 
2)7t), where k and / are integers. The matrix element 7 
for this tunneling process is very small, so that we are 
allowed to consider only the second-order processes with 
respect to 7. Thus, these lowest-order tunneling processes 
correspond to the change in ip p with m = ±1. Within this 
lowest-oredr approximation, Z± is simplified to 



Z± oc 1 ± 2cosx • Z 1 /Z , 



(46) 



where 

Zn 



I Vp p f 

J A(p p =2irm 



exp (-S 2 



V8V9 



'-'off 



S B ) . (47) 
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To calculate Z\, we first determine the stationary path 
of S st = S z + S A + S B under the condition of tp p (P) = 
f p (0) + 2ir, and then incorporate the influence of S^g, 
which plays the role of a dissipative environment, by in- 
tegrating out the low-energy fluctuations around the sta- 
tionary path. We introduce the new variables: 



Or 



9r = 



(m L (p p - MO) 



(48) 
(49) 



m L + M 

where tul = L/(2irv p K p ). Note that 8 r (Or) is the rela- 
tive (center-of-mass) coordinate with respect to ip p and 
—9. Then, S st becomes 

S st = J drl — (d T 9 r y + —{d T 9 R y + — (d T 9y 



4V . 6 

cos r cos - 

na 2 



(50) 



where M r = m L M/(m L + M) and M R = m L + M. 
The stationary path is determind under the boundary 
condition of 6 R (/3) = R (0) + 2TTm L /(m L + M) and 
r (/3) = r (0) + 2n. Noting that m L > M, we find that 
Or ~ <p p - Thus, we approximately obtain the stationary 
path as 



2tt 



flst St 

9r~<P p - jr. 



1 



0?= 2 / (r-ri) + -/(r-r 2 ) 

9 st - I(t - n) - I(t - r 2 ), 



(51) 

(52) 
(53) 



where I(r — Tj) describes one instanton at r = Tj and 
satisfies I(—oo) = and I(oo) = 2tt. Then, we find that 

^ = I/( t - Ti ) + 1/(t-t 2 )-^t. (54) 

We calculate Z\ taking account of the low-energy fluc- 
tuations around the stationary path. Due to the barrier 
potential, 9 and ip p + are strongly pinned. Thus, we can 
neglect the fluctuations of 9 around 9 st . In contrast, (p p 
and 9 can fluctuate under the condition that Lp p + 9 is 
strongly pinned. This means that they are expressed as 



<p p = <f$ + <j>, 
= 9 st - A 



(55) 
(56) 



where <fr represents the low-energy fluctuations. We ap- 
proximate that I(t) = 27t?9(t), where i9(r) is the step 
function. The Fourier transform of 9 st and that of 9 st are 
expressed as 



— 1 7T 17T • 

9 st (iu) = -e 1WTl + -e IWT2 , 



9 st (oj) 



i2n 



i2tt 



(57) 
(58) 



Substituting eqs. (57) and (58) into eq. (47), we find that 

|i =7 2 [ P dndT 2 ( V<j)e~ s ™, (59) 
Zq Jo J 
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irL 



a z^ 



Llu 2 



v p K p (3 f}^ 27T Vp K, 



-<f>(u)<j)(-(j) 



+ i^E^ 1 H^H^ t (-)- (60) 



Integrating out 0, we obtain 

Zy 
Z 



7 2 exp 



— — T) / dndTue-*", (61) 
k p v p J Jo 



where 

Sir 



y L 

Z^ _l_ ,.,2 



47T 

27T x-- / 1 



+ 



x (1 -cosw(ti -t 2 )). (62) 

We substitute eqs. (38) and (39) into eq. (62). Noting 
that L ^> £, we obtain 

- 4-7T 1 , / Lui 

V2^ 



pf A>w>0 



+ E -coscch f — ) (l - COSUJ ( T i - T 2)), 
K p p f^ n oj i 



u)>0 



(63) 



where we have introduced the high-frequency cutoff A in 
the first term. 

We first consider the high-temperature regime of 
Vf/L < T < A. Noting that 



^2 ex P ( ~P 



M = -oo 



M cx 1 + 2 exp 



TlKn-L , 



(64) 



£ exp(-^M 2 ).l-2exp(- 



M = 
(M=odd) 



(65) 



the partition function is expressed as 



Z cx 2 + 8 cos x • 7 exp — ir 







L K„L 



v p K p 



T 



L 



x / dndr 2 e 



"Sin 



(66) 



In the high-tempcrature regime, S^ns is obtained as 



K n 



■In 



2irT J 



(67) 



J. Phys. Soc. Jpn. 

Substitution of eq. (67) into eq. (66) yields 

'2ttT n 



Full Paper 



Yositake Takane 



Z oc 2 + 8 cos x ■ 7 /3 



A 



x exp — 7r 



v p K p 



+ 



K„L 



T 



(68) 



From eq. (26), we finally obtain 



j{x) 



x exp — 7r 



L 



K„L 



VpKp 



T)sm X . (69) 



Next, we turn to the low-temperature regime of T <C 
vf/L. In this regime, we are allowed to retain only the 
term with M = in eq. (27). Thus, the partition function 
is given by 



Jo 



Z oc 1 + 2 cos x ■ 7 / dndT2 e 
'o 

To evaluate Si ns , we approximate as 

47T 



(70) 



1 , /LaA 2 , /2cT AL\ ,„ , 
- tanh — rj— In , 71 

LO \2VpJ K p \ TT v p J V ' 



where 7 w 0.5772. The second term in S- ms is obtained 
as 



4-7T 



UJ>0 



— coscch ( ) (1 — cosw (n — r 2 )) 



A^ 



■In 



cosh(^(r 1 -r 2 ))]. (72) 



Using eqs. (71) and (72), we obtain 



Jo 



dridr2 e 



'A \Al) 



(73) 



where c is a numerical constant of order of unity. Substi- 
tution of eq. (73) into eq. (70) yields 



(74) 



smx- 



(75) 



We finally obtain 

This indicates that the critical current in the low- 
temperature regime behaves like j c oc (l/L) 2K p _1 . 
This result is not consistent with Fazio et al.'s result, 
j™° oc (l/L) K p + K '~ 1 . We discuss the reason of this 
inconsistency in the final section. 

4. Josephson Current in the Weak-Coupling 
Model 

In this section we calculate the dc Josephson current 
in a TL liquid of length L with open boundaries, which 
is weakly coupled with the left (right) superconductor at 
x = [x = L) through a tunnel junction (see Fig. 1(a)). 
We show that the resulting expression of the Josephson 
current is essentially equivalent to that in the interface- 
barrier model studied in §3. 



The bosonized description of a TL liquid with open 
boundaries is presented by Fabrizio and Gogolin. 29 In 
terms of the boson operators a q , a q , b q and b q , and the 
winding-number operators N and M, the Hamiltonian is 
given by 



H = 



AL 



K p K n 



-M' 



where q — nn/L (n = 1,2,3, 



{v p qa\a q + v a qb\b q ) , 
(76) 

• • ) . The winding- number 



operators satisfy the condition of N + M = even. The 
phase fields are given by 



e+(x) 


ttN 


+ 6+(x), 


(77) 


e_(x) 


= V P + t 


?_(*), 


(78) 




ttM 
= — X 


+ 4>+{x), 


(79) 






1>-{x), 


(80) 



where 6± and <f>± are the nonzero-mode components, and 
ip p and ip a are the zero-mode operators which satisfy 
[N, ipp] = [M, if a ] = 2i and [N, ip a ] = [M, <p p } = 0. The 
nonzero-mode components are given by 



0+(x) 
4>+{x) 



q>0 



9>0 



IK. 



^E 

q>0 



P q>0 



qL 
2T 



e- aq / 2 sm(qx) (a\ + a q ) , (81) 
c- aq ' 2 cos(qx) (oj - a q ) , (82) 
e~ a ^ 2 sm(qx) (b\ + b q ) , (83) 



IK 



'JLe- a ^ 2 cos(qx)(bl-b q ). (84) 



It is clear that sin(gx) = when x = and L, so that 
0+(x) and 4>+(x) are pinned at the ends of the ID system. 
That is, 6>+ (0) = 9+(L) = and 0+(O) = <j>+{L) = 0. 
This means that fluctuations in both charge and spin are 
strongly suppressed in the vicinity of both the ends due 
to the open boundary condition. This result plays an im- 
portant role in our argument. In terms of the phase fields, 
the right-moving component ipi s and the left-moving 
component ip 2s of the electron field are expressed as 29 

i>u(x) = -J= exp (ik F x + ^{6+(.x) + 0_(x) 



V2 



net 



V2 



ira 



■. exp i 



s$+(z) + s$- (85) 



ik F x + -{ - e + (x) + e-(x) 



s<&+(x) + s$_(x)} 



(86) 



We neglected Majorana fermions in eqs. (85) and (86), 
which do not play an essential role in our argument. 

Now we introduce the Hamiltonian which describes the 
superconductors and their coupling with the TL liquid. 
We assume that our STLLS system is symmetric with 
respect to the NS contacts. That is, the left contact at 
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x = is equivalent to the right contact at x = L. Thus, 
we present here only the part of the Hamiltonian related 
with the left superconductor. The Hamiltonian for the 
left superconductor is given by 



k,s 

E { AciXlc kA-K- + Ae- i ^c L _ fc ,_c Lfc , + ) , 
k 

(87) 

where e k denotes the single electron spectrum and CLfc, s 
( c Lfc s) i s the fermion anihilation (creation) operator. In 
terms of the Bogoliubov transformation 

CLfc,+ = u k d-L k ,+ + Vke 1Xl dl_ k _, 



CL-k- = Ukd-L-k- - v kQ 1Xl dl k + , 



with 



£k_ 

Ek 



the Hamiltonian is diagonalized as 



H s - y^E k d Lk s d LktS , 



(88) 
(89) 

(90) 
(91) 

(92) 



k.s 



where E k = \A| + A 2 . The weak coupling between the 
TL liquid and the left superconductor is described by the 
tunneling Hamiltonian 

HT = 7f E ( c ks f dxt Lk (x)Mx) + h.c\ , (93) 

k,s ^ '' ' 

where V is the volume of the superconductor. We assume 
that the tunneling-matrix clement t^ k (x) has nonzero 
values only in the vicinity of x = and vanishes when 
otT < x 7 where «t represents the length scale of the or- 
der of a few lattice spacings. The Hamiltonian for the 
right superconductor is simply given by replacing L with 
R in cqs. (92) and (93). Note that t-Rk(x) has nonzero 
values only in the vicinity of x = L and vanishes when 
x < L — cut- 

In order to calculate the partition function, it is con- 
venient to derive the effective action 5 r describing the 
coupling with the TL liquid and the superconductors. 
The action S r is obtained by integrating out the elec- 
tron fields in the superconductors, 

/ U Vd ks Vd ^sVdl k s Vd RKs cxp (s s - S T ) 

k,s 

oc cxp (S r ) , (94) 
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/ dT E [ d Lk,s d r d ^k,s + E k dl k s d Lk , s 
J ° k,, 

+ d Rk,s d T d Rk,s + £fc<4fc,s rfRfe ' s ) ' ( 95 ) 

S T = f Q dr^= J2 (cL, 5 / dxt hk (x)Mx) 

+ c R k ,sJ dxt Rk (x)ip s (x) +h.c.y (96) 

The derivation of S r is outlined in Appendix C. The 
result is 

S r = / dndnQin - r 2 ) 
Jo 

x (e- ixi {^i_(0,ri)^2 + (0,75) 

+ 1p2- (0,Tl)lpl+ (0,T 2 )} +h.C. 
+ C- i ^{^ 1 _(L,T 1 )V 2+ (i,T 2 ) 

+ ^ 2 _(L, r 2 )}) + h.c.) , (97) 

where Q(t) is defined in Appendix C. We substitute 
eqs. (85) and (86) into eq. (97). Noting that 0+(O) = 
0+(L) = and 0+(O) = <j>+(L) = 0, we obtain 

i r 

S r = / dr 1 dr 2 Q(r 1 - r 2 ) 

x e - ixi expi(e_(0,ri) + e_(0,r 2 ) 

-$_(0,Tl)+$_(0,T 2 )) +C.C. 

+ (-l) M e-^ exp i (e_ (L, ri) + 0_(L, r 2 ) 



$_(L,n) + $_(L,t 2 ) +c.c 



(98) 



The kernel Q(r) vanishes when |r| » A -1 as shown in 
Appendix C, while the characteristic time scale of 6+ 
and $_ is much longer than A -1 . Thus, we can approx- 
imate in eq. (98) that 



Q( T i - T j) = r<5 ( r i - T J")> 



(99) 



where r = J AtQ(t). Thus, S T is simplified to 

S T = -— j\r\U-^ +i@ -^+ C A 
naj [\ J 

+ (-l) M (e- i ^+ ie -^)+c.c.) 



(100) 



The sign of the second term depends on M, so we express 
S r with even M and that of odd M as S\ and SE, 
respectively. 

Note that S± does not contain <p a , so that M is a 
conserved quantity. By contrast, N is no longer conserved 
in the presence of S± . The partition function is expressed 
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as Using this decomposition, we obtain 



M 

(M =cvcn) 



(M = odd) 



with 



E GXP (^^tr) Z+ = E e "^ m2 / [ T\ValVa q 

+ I V, (!»!) x-P(-^-^), (.08) 

v 7 u m= — oo 

x e " : ^ £m2 e im ^±+ i2 ¥ £!l ( Tl - T2 ) 
x /"x>0„ / [] PaJPa, exp (-S z - S NZ ) 

A J J q>0 

(102) x < c lx cxpi^ p (ri) - (^ p (r 2 ) +0_(O,ri) 



z±= f v ^p I n pa 9 pa 9 



x exp [' dT ^ d ^ P - S 2 



where ( + = and = 7r, and 



/3 



r - 0_(L,t 2 ) N ) +c.c. I. (109) 



dT wH^< (103) 



j p^p 



To evaluate Zi±, it is convenient to introduce the corre- 
S NZ = J'" dr ^ (aj0 T a, + Upgaja,) . (104) lation function defined by 

-t 2 ) =^ expi^p(n) - ^ p (r 2 ) 



<J>0 

We introduced £± to maintain the constraint N + M 

even. 33 Since it is assumed that the coupling between the _|_ 0_(O, n) — 9-(L,t 2 ) ) ), (110) 

TL liquid and the superconductors is weak, we calculate 7 ' 

the partition function by treating S± as a perturbation. where 

Within the lowest-order approximation with respect to j T>ip p jY\ >0 T>a^T>a q ■ ■ -cxp (— S z — 5 NZ ) 

S£, we find that Z ± oc 1 + Z 1± /Z 0± , where <- ) = j^^^^Z^T^i) ' 

(111) 

The correlation function is obtained as 



Z °± = E / P ^ / ~[[ValVa q 

m =-oo jA Vp=27rm J q>Q 

/ nL 47T ^-^ 1 / 

x exp (im C± - S z - S NZ ) , (105) « M = CX P ( - ^— g - ^ E 



Zl±= E / A , ^/ll^4(^) 2 --^-Vlcosechf^Vl-coscr)Y (112) 

x cxp (imC± — S z — S NZ ) . (106) where the high-frequency cutoff A is introduced for the 
Using eq. (100), we find that second term in the exponent. By using O, we can express 

' Zi±/Z ± as 

Z 1± =± V e im <± [ V Vp [ T[Vap>a q Z 1± 2 [f> 

m ^ JA Vp =2, m J ^ — =±2 7 cosx^ dndr 2 

x exp (-S z - S NZ ) ^£ m2 imC±+i ^ (Ti _ T2) 

x — / dridr 2 E m e ^ 

^ 77(1 ; (113) 

where we set 7 = T/ (7ra). 

We first calculate the Josephson current in the high- 
(107) temperature regime of v p /L <C T <C A. In this regime, 
we are allowed to retain only the term with m = 0, and 
where \ = X2 — Xi an d we neglected terms which obtain 
are independent of \. To proceed, we decompose (p p as ^ 

^p( t ) = Vp + 27rmr//3 + ^p(T), where ^ p (/3) - <^ p (0) = 0. £i± = ±2 7 2 cos x / dndr 2 (n - r 2 ) , (114) 

^o± Jo 



x ^e ix expi(e_(0,n) - 9_(L,t 2 )) +c.c.^. 



where 
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The partition function is 



Z oc 1 + 2e" 



+ 1 - 2e" 



= 2 + 8cosx • 7 2 /? 2 



2ttT 



Zi+ 
Zq+ 



x exp — 7T 



+ 



K„L 



(116) 



We have used eqs. (64) and (65) in deriving eq. (116). 
The Josephson current is obtained as 

, \ „ ^ainZ 
j(x) = -2eT^- 



x exp — 7r 



+ _2_]T sin X . (H7) 



This result is equivalent to that obtained in §3 in the 
high-temperature regime. 

Next we consider the low-temperature regime of T <C 
Vp/L, where the terms with M — dominantly con- 
tribute to Z. We then find that Z cx 1 + Zi + /Zq + . Using 
Poisson's summation theorem, we rewrite the summation 
over m in Z\ + , 



E 



m— — oo 



4ttK p L 



(118) 

This equation indicates that we are allowed to retain only 
the term with q = 0. Then, we obtain 



Zo+ 



„ p 2 

= 2cosx-7 2 /3 2 / dridT 2 e _1 W (Tl_T2) 



with 



SI (r) w c' 



xfi(ri-7u), (119) 

'(ft)*-*®-)"*. (120) 

where c' is a numerical constant of order of unity. The 
partition function is obtained as 

Zcxl + 2cos X . C7 2 |(-^)* _1 , (121) 

where c is a numerical constant of order of unity. We 
finally obtain 

„ ^dlnZ 
j( X ) = -2eT— 

= 4 ^(i) 2 (lz) 2( "" l)si - ™ 

The result is also equivalent to that obtained in §3 in the 
low-temperature regime. 



5. Discussion 

We studied the dc Josephson effect in a TL liquid 
on the basis of the two different models. We derived 
a bosonized description of a superconductor-TL-liquid- 
superconductor (STLLS) system with perfect NS inter- 
faces, and calculate the Josephson current by introducing 
a strong barrier at each NS interface. We next exam- 
ined the Josephson current through a TL liquid with 
open boundaries, where the TL liquid is weakly cou- 
pled at its left (right) end with the left (right) super- 
conductor through a tunnel junction. It is shown that 
the two models provide us the same expression of the 
Josephson current in both the high-temperature regime 
of v-p/L <C T <C A and the low-temperature regime of 
T <C v-p/L. We found that the Josephson current in the 
high-temperature regime is given by 



x exp — 7r 



L 



K„L 



VpK p 



T | Mil \. 

(123) 



while 



in the low-temperature regime. 

We treated the strong-barrier limit of the interface- 
barrier model in §3. Let us briefly consider the opposite 
weak-barrier case. For simplicity, both the barriers are 
assumed to have the same strength. We employ the no- 
tations presented in §3. The action describing the weak 
barrier potential is 



S a = 



4F 



dr cos((p f 



i cos 



2' 



(125) 



with Vo/(aA) <C 1. We have used the fact that the spin 
degrees of freedom are frozen in the vicinity of the NS 
interfaces. The potential becomes large after renormal- 
ization. We apply a perturbative renormalization-group 
argument to obtain a scaling equation for the barrier 
potential. 6-10 If the band width A is reduced to fi, the 
potential becomes 

V ^) = V (A)^) Kp ~\ (126) 

as long as v p /L <C n < A. This means that the potential 
is renormalizcd to be large in the repulsive interaction 
cases of K p < 1. Thus, the system flows to the strong- 
barrier limit, where our analysis in §3 is justified. Since 
the coupling between a TL liquid and superconductors 
is reduced by the strong barriers in this limit, we expect 
that the weak-coupling model also describes correct low- 
energy physics. Our results given in §3 and §4 support 
this reasoning. 

To elucidate characteristic features of the interface- 
barrier model, it is instructive to compare our problem 
with the double-barrier problem in a TL liquid. Kane 
and Fisher, 7 and subsequently Furusaki and Nagaosa, 9 
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studied transport properties of a TL liquid in the pres- 
ence of a double-barrier structure, which consists of two 
5-function barriers at x = and L. This system is sim- 
ilar to our Joscphson system except for the absence of 
superconductors. The Hamiltonian of the double-barrier 
problem is given by H — H 1B + H B , where 




(127) 



H B = - 



4^o 



TXOL 



cos t) cos - COS COS — 
2 2 



+ sin 9 sin - sin (f> sin — 



(128) 



We assumed that k-pL = (mod ir). The variables 9, 9, cf) 
and 4> in eq. (128) are_defined as 9 = (9+(L) + 0+(O))/2, 
9 = 9+(L) - 9+(0), 4> = (ML) + 0+(O))/2 and = 
4>+{L) — 4>+(0), respectively. Note that 6L and arc 
not contained in H B . Thus, the partition function is 



where 
S 1D 



Z = j V9+ j V0+ cxp (-S 1B - S B ) 
= dr 

JO J-oo 



(129) 



4nK 



+ 



89+ 
dx 



1 



09+ 



f / dr 

JO J-i 




S B = 



f 

Jo 



AtH b (t) 



(130) 



(131) 



We can simplify the partition function by using the ef- 
fective action for 9,9,6 and 6. The result is 9 



where 



Z = j V9V9 J VSV^cxp - S B ) , (132) 

6(uj)6(-(j) 



1 



9 1D = 1 V - 

off o~ is n i 



2\u\ 



1 



+ 



+ 



+ 



+ exp(-L\u\/v p ) 
2U 



8irK p f3 
1 

2nK a f3 
1 



St 



ex${-L\u}\/v p ) 
gM 

+ exp(-L|w|/tv) 
gM 

- exp(-L|w|/tv) 



0(w)(9(-w) 
w)t/>(-w) 

(133) 



We can clarify the influence of superconductors on a TL 
liquid by comparing the effective action S^g given in 
cq. (37) with that of the double-barrier problem 
given in eq. (133). As noted in §3, we see from eqs. (40) 
and (41) that J~ 1 (o;) and J a r 1 (w) arc of the order of A 
in the limit of uo — > 0. This indicates that (f> and </> have 
a mass gap of the order of A in our problem, while such 
large gap does not appear in the double-barrier prob- 
lem. We can understand the feature by noting that only 
a spin-singlet pair of two electrons can transfer the NS 
interfaces due to Andreev reflection. The conservation 
of the total spin through the Andreev reflection process 
leads to suppression of the spin fluctuations near the in- 
terfaces. The mass gap of the order of A reflects it. The 
presence of superconductors also affects the charge de- 
grees of freedom. We also see that J~ 1 (w) and J^ 1 ^) 



are of the order of v p /L in the limit of i 



0. Thus, 9 



and 9 have a small mass gap of the order of vp/L. This 
gap is simply induced by a finite-size effect. Let us fo- 
cus our attention on the regime of v p /L <C |w|, where 
the finite-size effect plays no role. The part of S^g de- 
scribing the dynamics of 9(uj) with v p /L <C \oj\ is given 

by 

'^ = 7^ Y. \u\§(u)0(-u>). (134) 



cNZr 

J cff 



2irK p ^ 

P ' v p /L<^\(v\ 



Similarly, the part of S^g is 
1 \ 



^eff 



2irK p /3 ^ 

pl d p /L<|w| 



2\uj\e{uo)6(-uo). (135) 



This indicates that the fluctuations in 9 are enhanced 
by a factor of two in our Joscphson system compared 
with the double-barrier system. The fluctuations of 9 arc 
also enhanced by a factor of two. This feature is also 
attributed to Andreev reflection. The reason is that the 
unit of charge transfer is not e but 2e in the Andreev 
reflection process. This doubling of the unit charge en- 
hances the charge fluctuations in the vicinity of the NS 
interfaces. 

We obtained that the critical current is j c oc 
(l/L) 2 ^ 1 - 1 when T < v F /L. As noted in §3, our 
result is not consistent with Fazio et al.'s result, 19,20 
jf HO oc (l/L) K p + K °- 1 . The reason of this inconsis- 
tency is now clear. Fazio et al. treated a TL liquid of 
infinite length, which is weakly coupled with the left 
(right) superconductor at x = (x = L), under the 
assumption that the potential induced by the coupling 
with the superconductors can be completely neglected. 
Thus, within their assumption, the dynamics of the phase 
fields in the ID system is completely unaffected by the 
coupling. However, even a very weak potential becomes 
large after renormalization due to the repulsive electron- 
electorn interactions in a TL liquid, 6 so that the TL liq- 
uid is effectively disconnected at x = and x = L. In 
such a situation, 9+ and </>+ are pinned at the discon- 
nected points, while fluctuations in and are en- 
hanced around there. From the above argument, we see 
that Fazio et al.'s result may apply to only restricted sit- 
uations, where the induced potential is completely negli- 
gible. If the potential becomes very large after the renor- 
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malization, Fazio et al.'s model must be replaced with 
our weak-coupling model, as suggested by Fabrizio and 
Gogolin. 29 This means that our result is much general 
compared with Fazio et al.'s result. 

Details of the difference between j c oc (l/L) 2K p _1 



and j c tHO oc (l/L) K p 



k- 1 +k„-\ 



are explained as follows. 



The dynamics of the phase fields is assumed to be unaf- 
fected by the coupling with the superconductors in refs. 
19 and 20, while, in our weak-coupling model with the 
open-boundary condition, we explicitly incorporate the 
strong pinning of and the enhanced low-energy fluc- 
tuations of 6L in the vicinity of the NS interfaces. The 
correlation exponent K a does not appear in the expres- 
sion of j c due to the pinning of <p + , while HO contains 
K„. The prefactor of K~ l in the exponent of j c is dou- 
bled compared with j™° due to the enhancement in the 
fluctuations of 

Appendix A: Bosonization of TL liquid sand- 
wiched between superconductors 

In this appendix, we adapt the bosonization method 
to a TL liquid sandwiched between superconductors. We 
first bosonize the ID electron system in the noninteract- 
ing limit, and then incorporate the influence of electron- 
electon interactions. Our bosonization precedure is based 
on the argument by Haldane. 17 

In terms of the eigenf unctions of the BdG equation, 
presented in eqs. (2) and (3), the electron-field operators 
defined in eq. (9) are expressed as 30 



1>i+(x) = E ( u i( x ) c i:+ - 9*q( x ) c l,-) , (A-l) 
<?>o 

1Mb) = E (/*(*)<**,+ ~ . ( A '2) 

<?>o 

i>l-( X ) = E ( U q( X ) d q- + 9q( x ) d q,+ ) > ( A '3) 
<?>0 

V>2-(*) = E {fMc*,- - «»4+) . ( A4 ) 

<j>0 

where c q ^ s and d qiS (c q s and d q s ) are the fermion ambula- 
tion (creation) operators and q = ir(n+l/2)/(L+£) with 
n = 0, 1, 2, • • • . The linearized quasiparticle dispersion is 

2 



e q ±=v F [q± 



X 



1 



X 



(A5) 



2(L + 0/ 2mV2(L + 0, 
Thus, the quasiparticles in the ID system are described 
by the Hamiltonian 

,2 



Hn = 



V F X 



4tt i + C 



q>0 

E '' |; ' 



g>0 



2(^ + 0, 

^ (c\_C q - +4,+ d 9,+) 



2(£ + 0. 



(A-6) 

Here, the first term represents the %-dependent compo- 
nent of the ground-state energy, arising from the second 
term in eq. (A-5). 



We define the density operator Pj S {x) (j = 1,2) as 
Pj S (x) — ^ s (x)ipj S (x). Its Fourier transform is given by 



Pjs(q) = / dxe lqx p js (x), 



(A-7) 



where g = 7r(n+l/2)/(£ + £) with n = 0,±1,±2, • • • . An 
important role in our bosonization procedure is played by 
the new density operators, 

Pa (q) = c^ 2 Pl+ (q) - c-^ /2 p 2 -(-9), (A-8) 

Pb (q) = e*«/ 2 pi-(«) - c-^/ 2 p 2+ (-<?)• (A-9) 
For q > 0, they are expressed as 

Pa{q) = E ( c gi+g>+ c 9i'+ c 9i>- c 9i+g,-) 



9l>0 



-i E c -*i+*,+<4,-. ( A - 10 ) 



9>9l>0 



Pa(-g) = E (4l, + C 9i+^+ + C 9i+9:- C L-) 
<?l >0 



1 E C 9i,- C -9i+9,+ ' (A' 11 ) 
<3><3l>0 



Pb 



(?) = E ( < 4+d.- d 9i.- +d 9i. + < 4i+9, + ) 



9l>0 



-i E 4, 1+ «,-4,+> (A12) 



9>9l>0 



= E (4,- rf ?i+9>- +£ W?.+4,+) 



<?i>0 



i E d </i>+ d -9i+<z>-- (A-13) 

<j><ji >0 



We can show that 



[Pa{q),Pa{-q')\ = 



Lq 



lPb(q),Pb(-q')} = ~—Sq, 
[p a (q),p b (-q')}=0. 



(A-14) 



(A15) 
(A-16) 



The 9 = components, which characterize an excess 
charge with respect to the ground state, are given by 



Pa (q = 0) = N a + ^-, 
Pb (q = 0) = N b + ^- 1 



(A17) 
(A-18) 



where \/ {2tt) results from the shift of the occupied-state 
distribution due to the phase difference \- We call N a and 
N b the winding-number operators, whose eigenvalues are 
0, ±1,±2, • • • . We will employ the zero-mode operators 
(fa and ip b , which satisfy [N a ,ip a ] = [N b ,ip b ] = i and 
[N a ,y b ] = [N b ,<p a ]=0. 

Let us consider the commutation relations between the 
density operators with q ^ and the field operators pre- 
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sented in eqs. (A-l)-(A-4). We find that 




[Pa(q),tpi+(x)} 


= -e q{ +i/ >ipi + (x), 


(A19) 


[p a (q),i>2-(x)] 


= e- iq( ~ x+ ^tp 2 -(x), 


(A-20) 


[pb(q),i>i-( x )] 


= -e^^+^Vi-W, 


(A21) 


[Pb(q),ip2+(x)] 




(A-22) 


[pa(q),i>i-(x)] 


= [Pa{q),i>2+{x)] = [p b (q),tp 1 + 


(x)\ 




= [ Pb (q),^2-(x)}=0. 


(A-23) 
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From eqs. (A-31) and (A-32), we naturally find that 

pl + {x)=pl-{x) = ^(N a + ^). (A-34) 
We detcrmin so as to ensure the fermion commutation 



^J, ± (j,f = 1,2). The result is 



Noting eqs. (A-19)-(A-23), we express the field operators 

32 



+ -(N a -N b ), (A-35) 



as' 



^i s (a;) = 
■02s (ar) = 



1 



%-{x) 



'2-Ka 
1 



-ikFX-\-i&2s {x) 



V2 



-KOL 



(A-24) 
(A-25) 



-+z(^£)M) 



where 



M»)=ff+ + i I E — 



<?>0 



( c -i 9 (^?/2) pa((z) 



+ - (/V a - jV fe ) , (A-36) 

-I 

-^(N a -N b ), (A-37) 



_ e i,(x+£/2) po( _ g) ^ (A . 26) 



9 2 Z +(^) 



z(* + £) 



q>0 



-I 



-(/V a -jV fe ), (A-38) 



-"9/2 



<?>o y 



02- (x) = 6. 



where <f a and f b are the zero-mode operators mentioned 
above. Strictly speaking, if we need to ensure the fermion 
e lq ^Pb(—q)j, (A-27) commutation relation between tpj± and VVt an d that 

between ipj± and V'j'^i Majorana fermions must be in- 
q(x+i/2) troduccd in eqs. (A-24) and (A-25). However, since the 
Majorana fermions do not play an essential role in our ar- 
gument, we neglect them in the following. For actual cal- 
+ c~ lq ( x+ ^ 2 ) pb(— q) J , (A-28) culations, it is convenient to use the new phase fields: 31 

G+(x) = \(e 1+ {x) + 9!-(x) - 9 2+ (x) - 9 2 -(xj), 

(A-39) 



q>0 y 



+ 0" 



+i/2) Pa(q) 
iq{x+i/2) p a {-q)). (A-29) 



Here, a is a positive infinitesimal and 9j s (j = 1,2) rep- 
resents the zero-mode component to be determined. The 
phase field 9j s is related with the density operator pj s , 

«•<*> = h dJ t- (A - 30) 

For example, we find that 

p 1+ (x) = pl + (x) + (c-" (x+€/2) P-(9) 

q>0 

+ e i,(x+£/2) po( _ g )) ) (A-31) 



Q-(x) = \(0 1+ {x) + 0i_(z) + 2+ (x) + 9 2 -{x)) 

(A 

= \(0 1+ (x) - 9^{x) - 9 2+ {x) + 9 2 -{x)), 

(■ 

$_(a;) = - 0i-(x) + 9 2+ {x) - 9 2 -{x) 



(A-40) 

>)■ 

(A-41) 

■>)• 

(A-42) 

To express these phase fields compactly, we define f p , 
f a , J and M as 



p 2 _(x) = pUx) + ^J2( eiq 

q>a 



+ e- i9(x+£/2) Pa(-9)), (A-32) 



fp = fa + f b , 
fa = fa- fb, 

J = N a + N b , 
M = N a - N b . 



(A-43) 
(A-44) 
(A-45) 
(A-46) 



where 



P%(x) 



1 89 



3* 



2tt dx 



(A-33) 



It is easy to show that [J, ip p ] — [M, f a ] = 2i and 
[J, fa] = [M,ip p ] = 0. Substituting eqs. (A-35)-(A-38) 
into the above expressions and using these operators, we 
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find that AE(Nb) and the x _ dependent part of the ground-state 

9 + M = v , + 9 + W, (A-47) entgy ("P/4»)(X' 2 /L), we obtain 

e_w - 1 ( J+ s) (,+§)+._<*), (a-48) "» = ^ + s) 2 + ("' + £) > <**» 

We have expressed the field operators and the Hamil- 

&+(x) = —M ( x + L + — J + 4>+(x), (A-49) tonian in terms of the density operators. In order to sim- 

L V 2/ plify the expressions, we introduce operators for q > 0: 



$-(x) =<p a + (l>-(x), (A50) 
where 

-"9/2 / 

q>a 

x (Pa(<?) - Pa(-q) + Pb(q) - Pb(-q)) , (A-51) 

-a?/2 



M&) = i^E — — cos « [ x + 2 



(x) = 7 2^— — smq\x + - 



a q = ^ 




(A-62) 


a\ = y 




(A-63) 


Pq = \ 




(A-64) 




(*.(«)" /%(«))■ 


(A-65) 



x (Pa(q) + Pa(—q) + Pb(q) + Pb(-qj) , (A-52) They satisfy the boson commutation relation. Using the 

boson operators, we rewrite the phase fields defined in 



<f, + {x) = 1 ^ sin<z L + eqs. (A-51)-(A-53) as 

m*o = i J2 \/^ aq/2 cos «(*+§) K - a ") 



q>0 

x (Pa(q) + Pa(-q) - Pb(q) - Pb(-q)) , (A-53) 

, --9/2 / £\ (A-66) 



q \x + - 



<J>0 ^ " - / N \ ^ / ^ 



(,:) 



(A-67) 



x(Pa(«)-po(-?)-p5(«)+p6(-g))- (A-54) 9>0 

Using the phase fields, we rewrite eqs. (A-24) and (A- 25) 

ip ls (x) = e ifcFX+i(e + (x)+e-(x)+ s $ + (x)+s$-(x)) ; g>0 V 9 L V V 

V27ra (A- 68) 

(A55) 

^ 2s(a ,) = _J_ e -ifc F x+i(-e + (x)+e_(x)- s $ + (x)+ 5 *_(x))_ M&) = i E V^Z° ^ C ° S 9 ( X + I) ^ ~ ^ ' 

(A-56) (A-69) 

We note that the density operators with q ^ and if Similarly, the Hamiltonian is expressed as 
satisfy following relations 



Ho = TT (( j + -) 2 + m2 )+E Wf<? («S + W 



[p a (q),H ]=-V F qp a (q), (A57) 4L 



9>0 



[p 6 (g),flo] = - W 6(g)- (A-58) ( A ' 7 °) 

Note that since J = N a + N b and M = N a - N b , then J 

We separate the nonzero-mode component H^ z from ff - must be even (odd) if M is even (odd). That is, J + M = 

From eqs. (A-57) and (A-58) we find that even. 

NZ irv-p We take account of electron-electron interactions in the 

= ~ ^ yP^'P^-V + Pb(q)Pb(-q)) ■ (A-59) foUowing We emp i oy a mo del interaction Hamiltonian 

9>0 iJi nt given as follows: 

Next we consider the zero-mode component Hq . The en- „l 

ergy cost due to the excess charge A a (> 0) is obtained H int = 2ttv f / dx~S^ (g2S s . s ' + g^s. -»') Pis(a;)/02s' ( x ) 

as s.s 

A£;(tv q ) = J2 v F { ^^P + ^To) + ™f £ dx g (^ S;S , + ^ a ,_ a o 

« T?^ 2 + TT^a. (A60) x (pi a (x)pi.»(x)+p2.(a;)p2.'(a;)). 

IL IL ^ A . 71 ) 

This expression also holds for negative N a . Similarly, 
we obtain the energy cost due to Nb. Adding AE(N a ), 
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In terms of p a (q) and pb{q), the interaction Hamiltonian Finally, we rewrite 6± and <f>± in terms of a q , a q , b q and 
is rewritten as 



H\ nt , — — 



TTV F g 2 



+ ^yp- (Pa(q)Pa(-q) + Pb{q)Pb{~q)) 



Trv F g 4 



L 



J2Pa(q)pb(-q)- 



(A-72) 



We decompose Hi n t into the nonzero-mode component 
Hj^ and the zero- mode component H^ t , 

H? n t = -(92+ 9 ' 2 )J2 V -¥(alal + a q a q ) 



q>0 



q>0 



+ (.94 + g'i) ^ v ?q a q a q + (.94 - g'i) ^ w F9/?J/?g, 
g>0 q>0 

(A-73) 



w z — _ 



ttvf 

IT 

TTVF 



(g2 + g'2-g4-g' 4 )(j+^y 



+ ^^2-g , 2 + gA-g' i )M 2 . 



The nonzero-mode part H NZ = H ^ z + Hj^ is diagonal- 
ized as 

^nz = ^ (y p qa\a q + v a qblb q ) , 

q>0 

where 

Vp(a) = vf\J (i + 34 ± g'if - (52 ± g' 2 f- 

Here, the new boson operators are given by 
a q = a q cosh A p — a q sinh A p , 
b q = (3 q cosh \ a + (3 q sinh A CT , 

with 

52 ± .9 2 



tanh(2A p((T) ) = - 



+ .94 ± g'i 

Z — rrZ 1 rrZ 



(A-79) 



The zero-mode part H = + i?i£ t is expressed as 

flZ =s(»^( J+ 3 2+ fe M2 )' (A - 80) 

where ^ P ( CT ) is the correlation exponent 



K 



/ l + (. 94 ±^)-(g 2 ±^) 
1 + (54 ± .94) + (.92 ± ,9 2 ) ' 
Consequently, the total Hamiltonian is 



(A-81) 



H = 



AL 



^ (v p qa\a q + v a qb\b q ) . 

q>0 



(A-82) 



b q as 



+ (a;) = i^/K~ P ^2 \\^J^ aql2 cosq ( x 



q>0 



qL 



-aq/2 



smq x 



P q>0 



(A-83) 

(4 + a ?) . 
(A-84) 



0+(.t) = \/J e ~ Q9/2sin,Z { X + I) ^ + ^ ' 



-e - Q<?/2 cos 



(A-85) 
9^ + |) (b\-b q ). 



(A-86) 



Appendix B: Derivation of 

Let us consider Z NZ defined as 

Z xz = / H Va\Va q Vb\Vb q 

q>0 



x cxp (-5 NZ [{aj, a„ &J, &,}] - S B [{ a t, o„ 6j, 6J] ) . 

(B-l) 

(A- 74) We can rewrite Z NZ as 



^NZ 



(A75) 






xs(§- 


(A76) 






x6^4> 


(A-77) 


xS(j>- 


(A-78) 






x exp ( 



Y[ Va\Va q Vb\Vb q / VeV9V<f>V^ 

q>0 J 



(*-(*+(£)- MO))) 



x J V\ p V\ p V\ a V\ a J Y[ Va\Va q Vb\Vb q 



<expi^-S NZ [{ala q ,blb q }] 

+ [Wq, a qi b \, b q}, 0, V \). A<t, A~ CT ] 



(B-2) 



where 

V = I dr 



J dr X p (r)(. 



\ P (T)[e(r)--(e + (L,T) + e + (o,T)) 



+ x p (r)(e(r)-(e + (L,T)-e + (o,T)) 
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+ A~ ct (t)U(t) 



.(L,t)+0+(O,t)) 



+ A <T (r)U(r)-(0 + (L,r)-0 + (O,r)) 



cosh 



cosl, ( ^ ) I ). (13 T) 



(B-3) 

Note that 0+(L) and 0+(O) in cqs. (B-2) and (B-3) are 
functions of {a ql a q } as shown in eq. (14). Similarly, 
<p+{L) and <j)+(0) are functions of {b q ,b q }. Integration 
over a g , a q , 6 g and 6 J yields 

/ V\ p V\ p V\ a V\„ ( Y\Va\V aq Db\Vb q 

xcxp(^-5 NZ [{at,a 9 ,6t,6J] 

x +YP [{a\, a q , b\, b q }, 0, 0, <j>, 4>, A~ p , X p , \ a , \ a 



V\ p V\V\ a V\ a 



x exp 



irK, 



23 ^2 J p( u ) x p( u ) x p(.- u ) 

-il^MV-w) 



2-kK 



5^Jp(w)Ap(w)Ap(-w) 

-iI^0»A p (-a;) 



7T^ 

2/3 



J (7 (w)A(w) CT A (7 (-w) 



2irK a 



where 
J p (uj) 



LU 

-il^WKi-^y (B-4) 
cosh ( ) 



+ cosh 



2^sinh(^) 



+ cosh 



+ 1 



2« P 
(B-5) 



J p (w) - 



1 



2uj sinh ( ^ 

(L - €)w 



cosh ( ] 



+ cosh 



coshj ^ ) (B-6) 



1 



2w sinh ^ 



cosh 



— cosh 



+ cosh 



fa 



-1 . (B-8) 



Integrating out A p , A p , A CT and A CT , we find that 
J V\ p V\ p V\ a V\ a J Yl Va\Va q Vb\Vb q 



q>0 



exp (-5 NZ + YP) 



cx exp (S^) , 



(B-9) 



where 



^eff — 



1 ^J-V^w-U,) 



+ 



8tt# p /3 
1 

2wK a 3 
1 



^^'(w)^)^). (B-10) 



5 B [0,0, 0,0]). 



We thus find that 

Z NZ cx y V9V8V4>V0 

xexp(-S e N ff z [0,0 

(B-ll) 

It is clear that the effective action given in eq. (B-10) 
satisfies eq. (36). 

Appendix C: Derivation of S r 

We assume that our system is symmetric with respect 
to the left contact at x — and the right contact at 
x = L. For simplicity, we treat only the coupling at the 
left contact between the TL liqid and the left supercon- 
ductor, so that we do not explicitly express the subscript 
L, indicating the left superconductor, in the following 
derivation. 

We first rewrite the tunneling Hamiltonian given in 
eq. (93) using the Bogoliubov transformation, 



H 1 



= —Y 



k,s 



i 



2lo sinh 



< cosh f — — 



+ (u k J dxt k (x)tp + (x) 

+ v k c ixi J dxt*_ k (x)i>l(x)j 
- k -( u k J dxt- k (x)ip-(x) 
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dxt_ k (x)ip_(x) 
+ u k J dxt* k (x)^ + (xj^dk,- 
+ ( _ v k c~ 1Xl j dxt k (x)ip+(x) 

+ u k J dxt*_ k {x)^l{x))d. 



(C-1) 



The action S T is obtained by integrating out the electron 
field in the left superconductor, 

cxp (S r ) oc f IJl>4 >a X>dfc, a exp (S s - S T ) , (C-2) 



where 



S s = ^(-ie + E k )(dl + (e)d k , + (e) 

+ rf t _ fe; _(e)rf_ fe ,_(e)), (C-3) 

S T = / driJ T (r). (C-4) 
Jo 

Here, e denotes the fermion Matsubara frequency. After 
the integration, we find that 



u k v k 



P^V ^ -ie + E k 



x ! dxi_ fc (x)^-(a;, -e) J dyt k (y)ip+(y,e) 

+ e^| ^(^(x.e) | dyt*_ k (y)^_(y,-e) 

- c~ 1Xl J dxt k (x)ip + (x, -e) f dyt-k{y)i>-{y,e) 

-e ixi | dxt_ fc (a#l(*,e) | dy4(j/)^(y,-e))- 

(C-5) 

In deriving eq. (C-5), we neglected irrelevant terms which 
do not depend on the macrosocpic phase Xi- Equation 
(C-5) is simplified to 

r 1 \ 1 x 2E k 

x (e~ 1X1 J dx J dyt-k(x)t k {y)i)-(x,-e)i) + (y,e) 

+ h.c.Y (C-6) 



We assume that the tunneling-matrix element satisfies 

(t k {x)t- k {y)) kF = f(x)5(x - y), (C-7) 

where (• • • )fe F denotes the average over k on the Fermi 
surface. Note that F(x) has nonzero values only in the 
vicinity of x = and vanishes when ar < x. If we carry 



out the integrations over x and y in eq. (C-6) after averag- 
ing over k on the Fermi surface, the cross terms between 
ipi±(x) and -0 2T (x) vanish due to the presence of a spa- 
tially oscillating factor of e ±l2fcFa: . After the summation 
over k, we find that 



TA 



S3 V VA^ 



a (^i-(0,-e)^2+(0,e) 



+ ^2-(0,-e)Vi+(0,e)) + h.c), (C-8) 

where T = nN(0) J dxf(x) (N(0): density of states at 
the Fermi level). In the imaginary-time representation, 
S r is rewritten as 

5 r = £ dndnQfn - r 2 ) (^_(0, r^j+fO, r 2 ) 

+ V>2-(0,Tl)Vl+(0,T 2 )) +h.c.V (C-9) 



where 



TA 



(C-10) 



Note that we are interested in the case of T <C A, where 
Q(t) is approximated as 

TA 



i r 

Q{r) = - \ 

71" J A 



dz- 



-\r\z 



(C-ll) 



/a Vz 2 - A 2 

This indicates that Q{t) decays exponentially with a 
characteristic time scale of the order of A -1 . 

Adding the term describing the coupling at the right 
contact, we finally obtain 

S r = / dndnQin - r 2 ) 

Jo 

x ^e- ixi (v>i-(0,Ti)V2+(0,r 2 ) 

+ V>2-(0,Tl)Vl+(0,T 2 )) +h.C. 

+ e" iX2 (Vi_(L,ri)^+(£,75) 

+ ^ 2 _( J L,r 1 )V-i+(£,T 2 )) + h.c). (C-12) 
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